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Wc study the decay of quantum nonlocality, identified by the violation of the Clauser- 
Horne-Shimony-Holt (CHSH) Bell inequality, for two noninteracting Josephson qubits 
subject to independent baths with broadband spectra typical of solid state nanodevices. 
The bath noise can be separated in an adiabatic (low-frequency) and in a quantum 
(high-frequency) part. We point out the qualitative different effects on quantum nonlocal 
correlations induced by adiabatic and quantum noise. A quantitaive analysis is performed 
for typical noise figures in Josephson systems. Finally we compare, for this system, the 
dynamics of nonlocal correlations and of entanglement. 
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1. Introduction 

The presence of quantum correlations in composite nanosystems is an essential re- 
source for quantum information processing Amongst the most relevant aspects 
of quantum correlations, of both fundament al a. nd applicative role, are: entangle- 
ment 121^ quantum discord ^ and nonlocality Studying these quantities for re- 
alistic quantum systems that are promising candidates for realizing a quantum 
computer is an interesting issue. Considerable development has been recently made 
towards the implementatio n of a solid-state quantum computer. In particular, su- 
perconducting hi gh-fi delity single qubit gates with coherence times of ~ 1/iS 
are now available l-^. Two-qubit logic gates have been proved in different labora- 
tories ^ and Bell states have been prepared up to 75% of fidelity In a circuit 
quantum electrodynamic framework, highly entangled two-qubit states with con- 
currence up to 94% have been also generated, allowing the first implementation of 
basic quantum algorithms with a superconducting quantum processoJ^. 
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In order to process quantum information, it is important to establish how long 
the relevant quantum correlations can be maintained in noisy nanocircuits. Solid 
state noise may indeed represent a serious limitation towards this goal. Josephson 
junction-based experimental setups are often influenced by broadband and struc- 
tured noise, whose typical power spectra show a 1// low-freq uency b ehavior up to 
some cut-off frequency followed by a white or ohmic behavior I 12 | 13 | •j^'j^q presence 
of slow components in the environment makes the decay of the coherent signal 
strongly dependent on the experimental protocol being used EHI]^ Measurements 
protocols requiring numerous repetitions are particularly sensitive to the unstable 
device calibration due to low-frequency fluctuations. The main effect is dcphasing 
due to defocusing of the measured signal. Incoherent energy exchanges between 
system and environment, leading to relaxation and decoherence, occur at typical 
operating frequencies (about 10 GHz) where noise is white or ohmic. 

Low- and high-frequency noise affects quite differently single-qubit gates: adia- 
batic (low-frequency) noise typically lea ds to a lgebraic decay, while quantum (high- 
frequency) noise to exponential behaviorl ^^ l ^^ i It is therefore important to evaluate 
the effect on quantum correlations time evolution of adiabatic noise and of its in- 
terplay with quantum noise. This analysis has been recently addressed in detail for 
entanglement LtSl. Here we extend this study to the aspect of nonlocality. 

Nonlocal correlations, i.e. correlations not reproducible in the framework of Bell 
inequality tests by any c lass ical local model, are particularly important for quantum 
cryptography purposes ^1^. In this paper, we consider two noninteracting qubits 
subject to independent environments with typical broadband spectra. We use the 
Clauser-Horn e-Shi mony-Holt (CHSH) inequality and the maximum of the related 
Bell function to identify nonlocal correlations for classes of entangled initial 
states currently obtainable in laboratory. In particular, we study the sensitivity of 
the time evolution of nonlocal correlations to the initial state purity and degree of 
entanglement. 



2. Model 

We consider a system composed by two identical superconducting qubits, namely 
A and B, each interacting with independent baths characterized by a broadband 
spectrum. The total Hamilton ian is -fftot = Ha + Hb , where each single-qubit 

Hamiltonian is given by (h = ijnnn 

Ha ~ HQ_a — ^aC^z^af^, -ffQ,a = — ' i?q /2, (1) 

where Hg a refers to the qubit a = A, B, ffa is the Pauli matrices vector, {Vial = ^a 
the qubit frequency splitting and ^a are collective environmental variables whose 
power spectra are 1/f at / S [7™, 7m] (low- frequency noise) and white or ohmic 
at frequencies of the order of the qubit splittings (high-frequency noise). According 
to a standard model, noise w ith 1/f spectrum can be originated by an ensamble 
of bistable fluctuators (BFs) The physical origin of the fluctuators depends 
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on the specific setup. For instanc e, charge ba sed devices are extremely sensitive 
to background charge fluctuations I 14 | 19 | 20 [21J ^Qigg a,t higher frequencies instead 
may either originate from quantum impurities ^1^, possibly of the same physical 
origin, or being due to the circuitry. 

Effects of the different parts of the spectrum can be treated in a multi-stage ap- 
proach introduced in Ref.^^. Effects of low- and high-frequency components of the 
noise are distinguished by decomposing — ^ ia{t) + ^f.a Stochastic variables 
S,a{t) describe low-frequency 1// noise, and can be treated in the adiabatic and lon- 
gitudinal approximation. High-frequency (a; ^ fi^) fluctuations ^/.q arc modeled 
by a Markovian bath mainly leading to spontaneous decay. Therefore, populations 
relax due to quantum noise (Ti-type times), which also leads to secular dephasing 
(T2 = 2Ti-type). Low- frequency noise provides a defocusing mechanism determin- 
ing further coherences decay. In the Hamiltonian of Eq. ([T|) both the operating point 
(angle 9a between z and 17^) and the qubit splitting fie arc tunable. In the follow- 
ing we will consider both qubit operating at the optim al working point, 9^ = 7i'/2, 
where partial reduction of defocusing is achievecPH^. In addition we will focus on 
the case of identical qubits Vl^ = ^■ 

The two-qubit density matrix elements will be evaluated in the computational 

basis B = {|1) = |11),|2) = |01),|3) = |10), |4) = |00)}, where for each qubit 

we have i?Q,a|G) = — §|0), i/g^all) = §|1). Each subsystem "qubit -l-environment" 

evolves independently so that, once known the single-qubit dynamic J^, the evolved 

two- qubit density matrix will be readily obtained by a procedure reported in 
RefsPUl 



3. Maximum of the Bell function 



In this section we report the expression of the maximum of the Bell function for the 
class of two-qubit states whose density matrix px , in the standard computational 
basis B, has a "X structure" , i.e. 



PX 



/pii pu\ 
P22 P23 
P23 P33 
\pli paJ 



(2) 



This class of states is general enough to include the most common two-qubit states, 
like Bell states (pure two-qub it maximally entangled states) and Werner states 
(mixture of Bell states) ^ j^ ^ ^^ i Such a X form density matrix may arise in a variety 
of physical situations'^. A further remarkable feature of X states is that, under 
various kinds of dynamics, their s tructu re is maintained during time evolutior l^^ l ^^ l. 

Using the Horodecki criteriorff^l^, the maximum of the Bell function can be 
expressed in terms of three functions ui, U2 and of the density matrix elements 
as B = 2^u\ayijyk{uj -\- U/t}, where j, fc = 1, 2, 3. The CHSH inequality reads like 
B < 2, so that no classical local models are admitted for states such that B is larger 
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than the elassical threshold 2. The three functions 

Ul = 4(|pi4| + |p23|)^, W2 = (Pll + P44 ~ P22 ~ Pssf, "3 = 4(|pi4| - \p23\f ■ (3) 

Being ui always larger than ua, the maximum of the Bell funetion for X states is 

B = max{Bi, B2}, Bi = 2Vui +^2, B2 = 2y/ui + U3. (4) 

This quantity has been already studied in dynamical contexts of independent 
qubits each coupled to a bosonic reservoir (cavity) with Markoviaii^SI and non- 
MarkoviarPSI characteristics. In the following, we shall investigate the maximum 
of the Bell function, B, for our system of two independent Josephson qubits each 
interacting with individual baths. 

4. Initial states 

We consider extended Werner-like (EWL) two-qubit initial statei^H 

p*-r|$)($| + (1-0/4/4, /)*=r|«')(*| + (1-0/4/4, (5) 

whose pure parts are the one/two-excitation Bell-like states |(E>) ~ a\01) + 6|10), 
l^*) = a|00) + &|11), where |ap + |5p = 1. The purity parameter r quantifies the 
mixedness and a sets the degree of entanglement of the initial state. The density 
matrix of EWL states, in the computational basis, is non-vanishing only along the 
diagonal and anti-diagonal (X form) and this structure is maintained at f > in 
the system we are considering. Using concurrencJ^ C to quantify entanglement, 
one can also notice that the initial entanglement is equal for both the EWL states 
of Eq. © and reads C*(0) = C*(0) = 2max{0, {\ab\ + l/4)r - 1/4}. Initial states 
are thus entangled for r > r* ~ [1 + 4|a6|)^^. 

Entangled states with purity y 0.87 and fidelity to ideal Bell states « 0.90 have 
been experimentally generated^. These states may be approximately described as 
EWL states with rexp ~ 0.91. 

5. Time behavior of the maximum of the Bell function 

In this section we shall analyze the dynamics of the maximum of the Bell function, 
B of Eq. (U) , initially considering the case when only adiabatic noise is present and 
in a second stage including the effect of quantum noise. We shall also compare the 
dynamics of B and of the concurrence C. 

5.1. Adiabatic noise 

The effect of low-frequency noise components is obtained by the Hamiltonian 
Eq. ([1} with ^ « ^(t) treated in the adiabatic and long itudin al approximation, 
for which single-qubit populations do not evolve in timjm. The main effect 
of low-frequency fluctuations is dephasing due to defocusing processes. The lead- 
ing order effect depends only on the noise variance, which can be estimated by 
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independent measurements of the amplitude of the 1// power spectrum on the 
uncoupled qubits, as follows S^^^{uj) = 7rS]^[ln(7jvf /7m) oj]^'^- By exploiting the 
single-qubit c oheren ces determined in this case^^, we can construct the two-qubit 
density matriji^SUSI time t. Choosing initial states of the form ^ B is obtained 
from Eq. (j4]) . 

We find that, under adiabatic noise, B{t) = i?ad(0 is the same for both the 
initial EWL states of Eq. ([5]). The explicit expressions of -Bad(0 and of the times 
when Sad — 2, that is when a Bell inequality "violation sudden death" (VSD) 
occurs, arc given by 

R m 9 l^ I ,ad ^ / 4|ab|V^ . 

B^,{t) = 2r^l + ^^^^, t^^^^-^^^—^-l. (6) 

For any r < 1, the Bell inequality is violated at a finite time. Instead, for pure initial 
entangled states (r = 1), Bad{t) approaches asymptotically the classical threshold. 
The evolution of Bad{t) is displayed in Fig. [1] for = 0.02, which is a typical 



(a) (b) 




Fig. 1. Maximum of the Bell function, B, a.t = tt/2 and E/H = 0.02 under adiabatic (low 
frequency) noise. Panel (a) B as a function of dimensionless time Qt and |ap (r = 0.9); Panel (b) 
B as a function of dimensionless time f2t and r (a = l/\/2). 



figure of 1// noise in single-qubit experiment j I I [ The dependence on the initial 
degree of entanglement, parameterized by |ap, for fixed r = 0.9, is reported in panel 
(a). Sensitivity to the initial purity r, for fixed a = 1/^/2, is shown in panel (b). We 
note that the maintenance of nonlocal correlations {B > 2) strongly depends on the 
purity of the initial state r, while the dependence on |ap is weaker and symmetric 
around |ap = 1/2. The "violation sudden death" of Bad{t) for r < 1 is also clearly 
visible. 
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5.2. Interplay of adiabatic and quantum noise 

We now analyze the interplay of adiabatic and quantum noise simultaneously af- 
fecting the two units. 

When quantum noise (high-frequency, / ~ fi) is included, it adds up to the 
defocusing channel leading to both extra exponential decay of the coherences and 
evolution of the populations. This is due to the fact that, in practical situations, 
decay rates are much less sensitive than phases to fluctuations of control parame- 
ters. Singlc-qubit populations are obtained by the Born-Markov master equation. 
They decay with a relaxation rate T^^ = Sf (Q) /2 to asymptotic values which de- 
pend in general on 51 and on temperature jLL^ . For typical values of ft 10^^ 
rad/s) and T 0.04 K), asymptotic values of the excited, |4), and of the ground 
state, |1), populations are practically zero and one, respectively. O n the o ther hand, 
the coherences acquire an additional exponential decayim r factoj™. Using the 
reported single-qubit de nsity m atrix elements for this caseP^, we construct the new 
two-qubit density matrijiPSCS] from which we determine B{t). The maximum of the 




0.6 0.7 0.8 0.9 1 

r 



Fig. 2. Dependence of the VSD time on the purity r (a = l/\/2) for initial state p*. The 
behavior for p* is also quantitatively similar. The blue dashed curve is nty'gj-,, red dot-dashed 
curve is for quantum noise, black curve is the result of adiabatic and quantum noise altogether. 
Noise parameters are S = 0.02r2, Sf(uj) = 2 X 10*^ s'^. In addition, U = 10^^ rad/s, 9 = -it/2, 
T = 0.04 K. The inset zooms the region where t* 1, The point .Pcxp corresponds to T'cxp 

^ 0.91 

where CivsD ~ 3350. 



Bell functions for the two initial EWL states are now formally nonequivalent, this 
is a qualitative difference with the adiabatic noise case. However, for the typical 
experimental parameters involved in such nanosystems, quantitatively they do not 
differ significantly (their difference being always < 10^"^ in the violation region). 
We are interested to the VSD times ivsD when 5 = 2. These times are plotted 
in Fig. [5] as a function of the purity r for fixed a = l/V^ (maximally entangled 
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pure part), with noise parameters having values retrieved in experiment al [ In 
particular, the white noise level is 5/ = 2 x lG^s~^. We distinguish the cases of only 
adiabatic noise, only quantum noise and their interplay. For r < 1, adiabatic noise 
suppresses nonlocal correlations on a much shorter time scale than quantum noise. 
For high purity levels, the VSD time due to quantum noise is instead shorter than 
the adiabatic VSD time (inset of Fig. [2]), which goes to infinity for pure states (see 
Eq. ([6])). However, a quantitative estimate of the amount of nonlocal correlations 
preserved before B = 2 indicates that, for typical amplitudes of 1// and white 
noise, adiabatic noise considerably reduces the amount of nonlocal correlations on 
a short time scale even for ?' — > 1 . A finite value of tysD is ensured for any initial 
state only because of quantum noise. 



5.3. Maximum of the Bell function versus concurrence 
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Fig. 3. Maximum of the Boll function B versus concurrence C for the two initial Bell states 
|<E>) (blue solid curve) and I't) (red dotted curve) in the presence of both adiabatic and quantum 
noise. The points on the curves labeled with i = 1, 2, . . . , 5 indicates the values of (C, B) at times 
10~^f!t = i. B for 1$) decays a little bit faster than that for Noise parameters as in Fig. [2] 



Here we compare the dynamics of the maximum of the Bell function B with the 
evolution of the concurrence C, so to establish their connection in such a solid-state 
system. We consider initial preparation in the two Bell states |$) = (|01) + 1 10))/-\/2 
and l^") = (|00) + |ll))/-\/2. The behavior for different initial values of a and r does 
not differ qualitatively. Remarkably, for the considered system there is a one-to- 
one correspondence between B and C during the dynamics for both initial states, 
as sho wn in F ig. [31 This property has already been observed in other physical 
context J^^l^ (atomic qubits in cavities). However, the common behavior for both 
initial states is not predictable a priori. In general, it may strongly depend on the 
physical system and on the specific initial state 
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The slightly different time dependence for the two Bell states is evidenced by 
the different distance along the curves of the dots pointing to times lO^'^fJt = i 
{i = 1, . . . , 5). In addition it is clearly visible the threshold value of C below which 
there is no violation anymore. Starting from the initial state |$) (|^)) we find that 
for C < 0.43 (C < 0.38) the maximum of the Bell function _B < 2, so that we 
cannot be sure of the presence of nonlocal correlations in this region. 

6. Conclusions 

In this paper we have investigated the time evolution of nonlocal correlations (non- 
locality), identified by the maximum of the Bell function B when it violates the 
CHSH-Bell inequality {B > 2), between two noninteracting Josephson qubits sub- 
ject to independent baths with broadband noise typical of solid state nanosystems. 
In particular, an adiabatic (low-frequency) noise and a quantum (high-frequency) 
noise can be distinguished. We have shown that, while adiabatic noise has the main 
effect on nonlocality decay, it is the quantum noise that induces a complete dis- 
appearance of quantum nonlocal correlation for any initial state (even for pure 
maximally entangled states). In particular, we have also reported the times when 
B = 2 (ivso), after which there is no more certainty of the presence of quantum 
nonlocal correlations. We have finally compared, for this system, the dynamics of 
nonlocal correlations and that of entanglement, quantified by the concurrence C. 
We have found that a one-to-one correspondence between B and C occurs in time, 
independently on the form of the initial two-qubit state. Moreover, we obtained 
thresholds values of C below which the Bell inequality is not violated {B < 2). 

The results presented in this paper provide new insights towards the possibility 
to exploit nonlocal quantum correlations for quantum information processing with 
superconducting nanocircuits. 
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